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Acoustic  propagation  loss  is  o  critical  factor  in  the  design,  performance 
and  effectiveness  of  all  submarine  and  surface-ship  sonar  systems.  The 
attenuation  term  in  the  sonar  equation  for  propagation  loss  is  token  here 
to  include  all  the  mechanisms  for  which  the  loss  is  proportional  to  range. 
Absorption  in  the  medium  is  generally  the  dominant  mechanism. 
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Figure  1:  Mediterranean  propagation  loss  data. 

Absorption  tends  to  increase  as  the  square  of  the  frequency  end  this 
limits  the  maximum  useable  frequency  for  any  given  range.  The  situation 
Is  Illustrated  by  the  results  of  a  recent  sound-channel  experiment  in  the 
Mediterranean  Sea  111.  The  curves  of  propagation  loss  in  dB//im  shown  in 
Figure  1  clearly  demonstrate  the  increasing  effects  of  absorption  with 
Increasing  signal  frequency. 

In  order  to  deal  quantitatively  with  limitations  imposed  by  attenuation, 
the  other  loss  mechanisms  can  not  be  neglected  Internal  scattering  and 
diffraction  loss  ere  important  in  some  of  the  experiments  to  be  discussed 
However,  they  are  incidental  to  the  main  purpose  of  the  report,  which  is  to 
develop  a  predictive  global  model  for  sound  absorption  in  sea  water  that 
Includes  the  dependence  on  environmental  factors. 
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Background 

The  1935  Novel  Research  Laboratory  report  by  Stephenson  (2)  appears  to 
be  the  first  to  recognize  that  sound  absorption  in  sea  water  is  very  much 
greater  than  in  fresh  water.  However,  it  was  not  until  after  WWII,  after 
the  anomaly  had  been  clearly  established  by  extensive  field  measurements, 
thot  efforts  were  undertaken  to  determine  the  cause. 

Laboratory  studies  in  the  MHz  range  hod  shown  that  pure  water  followed 
the  classical  low  of  absorption  proportional  to  frequency-squared,  but  the 
magnitude  was  four  times  greater  than  thot  predicted  by  shear  viscosity 
clone.  This  onomoly  wos  eventuolly  ottributed  to  the  odditionol  effects  of 
bulk  viscosity.  When  extropoloted  to  the  10-100  kHz  range,  the  predicted 
coefficients  were  still  lower  then  the  sea  data  by  on  order  of  mognitude. 
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Figure  2:  Liebermonn’s  experiment. 

In  1948,  Liebermonn  13)  reported  experiments  comparing  fresh  ond  seo- 
water  absorption  in  the  band  10  kHz  to  l  MHz.  Deto  oblaineo  in  a  fresh¬ 
water  reservoir  agreed  with  high-frequency  laboratory  measurements. 
Below  100  kHz  the  sea  data  were  also  in  agreement  with  reported  ocean 
results;  however,  they  fell  off  at  higher  frequencies,  opproaching  the 
fresh-woter  volues  neor  1  MHz  os  shown  in  Figure  2.  Liebermonn  proposed 
ionic  relaxation  of  NaCl  as  a  candidate  mechanism. 
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Reloxatrunel  absorption  can  occur  in  aqueous  solutions  as  the  result  of 
pressure-dependent  chemical  equilibria.  Density  of  the  medium  changes 
according  to  LeChatelier’s  principle;  i.e.  if  cne  species  is  smaller  than  the 
other,  increasing  pressure  will  drive  the  reaction  in  the  direction  of  the 
smaller  species  and  the  compressibility  is  thereby  a'tered.  The  rote  of  the 
process  is  governed  oy  the  time  constant  of  the  reaction.  Therefore,  the 
sound  speed  c  becomes  complex  and  acoustic  waves  decay  with  r  ange 
Better  insight  is  gained  from  the  complex  wavenumber  k=2rtf/c.  Since 
dispersion  is  very  small  («1  m/s),  it  con  be  neglected  and  the  wavenumber 
can  be  approximated  as  ks  2Ttf/c0  +\A  where  c0  is  constant.  Then  we  have 

A=orfr  f2  /{i2  ♦  fr2)  where  A  is  the  loss  per  unit  distance  and  f  and  p^  are 
the  relaxation-frequency  end  amplitude  parameters. 

If  -an  shown  that  A„=  \  fr  *  k*  [Zl  AV?  pc  /2R(T+273)  where  k>  is  the 

forward  rate  constant,  W,  is  the  concentrafioi  of  the  species  Z,  AV  is  the 
volume  change  of  the  equilibrium  involved,  R  is  the  "gas  constant”  andT  is 
the  temperature  in  °C.  The  equilibrium  constant  for  the  reaction  is  given 
by  K=kV  ky  where  k<  is  the  rove-se  rate  constant. 

Various  types  of  equilibria  may  be  involved  and  each  one  has  a  specific 
formula  for  relaxation  frequency  and  volume  change  The  examples  listed 
below  are  typical  and  can  easily  be  extended  to  other  cases. 

The  simplest  type  is  the  single-step  equilibrium: 

Z  44  T  1rz(k>  *k*  )/2tj  K  =|Zl/iZ}. 

In  ionic  solutions  like  sea  water,  equilibria  can  be  of  the  type: 

Z«X**Y“  uWMYIWji  K={Z)/(X*)in 

Generally,  the  single-step  equilibria  are  diffusion-controlled  and  too 
fast  to  cause  sound  absorption  in  the  sonar  frequency  range.  Slower  ionic 
relaxations  in  sea  water  can  usually  be  modeled  as  two-step  equilibria, 
which  have  two  distinct  relaxation  frequencies.  The  second  step  is  rapid 
and  the  associated  high-frequency  relaxation  does  not  contribute  to  the 
absorption  at  sonar  frequencies  and  the  slow  first  step  controls  the 
relaxation  frequency  in  question.  This  type  can  be  represented  as: 

Z  o  Z*  X*  ♦  Y" 

Then  ff  =  (k*,*  k*,q/(  1  ♦q))/2n  and  AV  =  AVf  *  AV,  /(l*q) 
where  q=  K2  (IX*J  ♦  (Y'J),  Kt=|Z)/i2)  and  K2=  tZ!/lX*llV) 

Another  type  that  will  be  encountered  is  the  exchange  equilibrium 
XY~  e  YX*  tr  =  (k’dXMY'D*  k'ClYHX'm/ZK  K=lXilY")/IX‘)lY) 

The  formula  for  A^  remains  the  same  for  all  equilibria  except  that  the 
concentration  factor  !Z}  becomes  IXHY‘1  for  the  exchange  type 
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Figure  3:  Resonator  apparatus. 


In  1949,  Leonard  et.  al.  [41  reported  initial  results  of  a  laboratory  study 
to  determine  what  equilibrium  is  involved  in  sea-water  absorption.  Their 
apparatus  consisted  of  a  spherical  resonator  as  illustrated  in  Figure  3.  The 
resonator  was  suspended  in  a  vacuum  chamber  by  wires  to  minimize  any 
residual  losses  caused  by  acoustic  radiation.  The  useable  frequency  range 
was  25-500  kHz  and  values  of  Q>105  were  realized  in  pure  water  at  lower 
frequencies.  - 

In  the  resonator  method,  the  transducer  can  serve  as  both  transmitter 
and  receiver.  The  sphere  is  excited  specific  resonant  frequencies.  Radi8l 
resonances  usually  give  the  lowest  residual  losses  and  they  can  usually  by 
located  by  calculation.  However,  because  the  resonator  is  not  perfectly 
spherical,  the  resonances  are  hardly  ever  simple.  The  resonant  frequencies 
shift  with  changes  in  sound  speed  caused  by  variation  of  temperature  or 
the  addition  of  chemicals.  Interference  effects  among  adjacent  modes  can 
result  in  rapid  changes  in  residual  loss  and  this  limits  the  accuracy  of  the 
measurements. 

Losses  are  determined  by  switching  the  transducer  from  the  transmit  to 
the  receive  mode  and  measuring  the  rate  of  decay  in  dB/s.  The  system  is 
first  calibrated  using  distilled  water.  The  excess  losses  of  the  various 
chemical  constituents  are  then  determined  by  subtracting  the  calibration 
values  from  the  measured  values.  The  losses  in  dB/km  are  calculated  bu 
dividing  by  the  sound  speed. 

The  experiments  of  Leonard  et.  al.  showed  no  effects  for  NaCl  alone.  The 
salt  responsible  for  the  excess  absorption  in  sea  water  was  determined  to 
be  magnesium  sulfate  (f1gS04). 


Figure  4:  Resonator  measurements  of  Wilson  8nd  Leonard. 

Results  of  a  more  detailed  study  of  magnesium  sulfate  by  Wilson  and 
Leonard  (51  appeared  in  1954  and  are  illustrated  in  Figure  4  The  earlier 
conclusions  were  confirmed  end  dependence  on  concentration  was  also 
investigated.  It  was  found  that  Interactive  effects  of  NaCI  actually  causes 
a  reduction  In  sea  water  absorption 

In  1953,  Kurtze  and  T8mm  I6J  also  reported  results  of  similar  resonator 
experiments  confirming  Leonard  s  findings.  Several  pcssible  mechanisms 
were  proposed;  however,  it  was  not  until  1962  that  a  definitive  analysis 
by  Eigen  and  Tamm  171  disclosed  the  details  of  the  multistep  equilibrium 
involved.  Two  or  more  steps  were  indicated  by  the  fact  that  the  relaxation 
frequency  Increases  with  concentration  up  to  a  limiting  value. 

The  essential  acoustic  properties  of  the  relaxation  can  be  modeled  in  a 
more  simplified  form  as  the  two-step  equilibrium: 

(MgS04),  <->  (MgS04)„  «  Mg2+  ♦  S042- 

The  slower  first  step  involves  only  a  structural  change  and  this  controls 
the  relaxation  frequency  of  the  absorption  process  in  question.  The  second 
relaxation  is  very  rapid  and  its  contribution  is  negligible  The  reduction 
effect  of  NaCI  con  be  attributed  to  the  coupling  effects  of  another  rapid 
equilibrium  involving  the  ion-pair  NaS04*. 
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In  1962,  Schulkin  ond  Marsh  (8)  introduced  the  first  procticol  formula 
for  sound  absorption  in  sea  water.  The  approximation  was  based  on  the 
laboratory  resonator  experiments  and  extensive  sea  measurements  in  the 
frequency  range  2-25  kHz. 

The  parameters  involved  in  the  S&M  formula  ere  temperature,  salinity 
end  pressure.  The  pressure  factor  of  magnesium  sulfate  was  assumed  to  be 
the  same  as  the  pure  water  value  obtained  by  laboratory  measurements  at 
high  pressures  and  is  given  by  (1-6.54x1  O^P)  where  P  is  the  ambient 
pressure  in  atmospheres.  An  alternative  approximation  is  where  D 
is  the  depth  in  km.  The  S&M  fonnula  can  then  be  wriiten: 

A={0.02  S  fr  f2/(  fr2  ♦  f2 )  ♦  0.03  f2  /fr}x  10'D/35  dB/km 

fr  =59x1 0™9  kHz  I 

where  fr  is  the  relaxation  frequency  of.  magnesium  sulfpt'e,  5  is  salinity  in 

parts  per  thousand  (ppt)  and  T  is  temperature  in  °C.  The  first  term  is  the 
magnesium  sulfate  component  and  the  second  is  the  pure  water  component. 
The  temperature  dependence  of  pure  water  hss  been  approximated  by  using 
the  relaxation  frequency  as  a  parameter,  which  makes  both  coefficients 
equal 

For  the  parameter  values  S=35,  T=4°C  and  D=0,  the  formula  becomes: 

A s  50  f2/(  (7 1)2  ♦  f2)  ♦  0.0043  f2  dB/km 

which  reduces  to  A=0.01  f2  dB/km  for  f«  fr~7 1  kHz. 

The  S&M  formula  was  adequate  for  the  sonar  problems'of  that  period; 
however,  attention  in  underwater  acoustics  subsequently  shifted  to  longer 
ranges  and  lower  frequencies  and  the  accuracy  of  extrapolating  beyond  the 
data  base  was  questioned  A  number  of  sea  experiment?  were  then  carried 
out  which  indicated  anomalous  behavior  at  much  lower  [frequencies 
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Figure  5:  Thorp's  data  and  2-component  model. 

In  1965,  Thorp  191  carried  out  a  sound-ch8nnel  propagation  experiment  In 
the  Bermuda-Eleuthera  area  over  a  500  km  track.  Measurements  were  made 
at  1200  m  near  the  sound-channel  axis  using  explosive  sources.  Results 
showed  a  clear  anomaly  below  3  kHz  as  seen  in  Figure  5.  Comparison  with 
other  earlier  data  confirmed  his  conclusion  that  the  values  were  an  order 
of  magnitude  greater  than  that  predicted  by  the  Schulkin-Marsh  formula.  In 
a  later  paper,  Thorp  and  Browning  [  10)  modeled  the  anomaly,  adding  a  I  kHz 
relaxation  component  to  the  absorption  formula.  The  formula  used  in  the 
data-fit  of  Figure  5  is  given  by: 

A=  0.009  f2  ♦  0. 1 1  fr  f2/(  fr2  +  f2 )  dB/km  (fr=  1  kHz). 

The  first  term  is  the  magnesium  sulfate  component  approximated  for 
frequencies  well  below  the  relaxation  frequency  71  kHz.  The  second  term 
is  the  anomalous  component.  The  pure  water  component  is  negligible  in 
this  frequency  range  and  has  been  omitted. 
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Figure  6:  Leroy’s  data  and  2-component  model. 

Another  propagation  experiment  was  carried  out  in  the  Provencal  Basin 
of  the  Mediterranean  Sea  by  Leroy  11 1).  Ir  ihis  case,  propagation  loss  was 
measured  over  a  single  refraction  path  and  not  in  the  sound  channel.  The 
range  was  therefore  limited  to  40  km. 

Leroy's  results  also  showed  a  similar  anomaly  of  roughly  the  same 
magnitude  as  Thorp  but  with  a  significantly  higher  relaxation  frequency. 

The  formula  used  by  Leroy  for  the  data-fit  is  given  by: 

A=  0.006  f2  ♦  0. 1 55  fr  f2/<  fr2  ♦  f2 )  dB/km  (fr=  1 .7  kHz) 

The  first  term  is  the  magnesium  sulfate  component  and  the  coefficient  is 
in  reasonable  agreement  with  the  S&M  formula  for  the  1 3°C  temperature. 
The  coefficient  of  the  anomaly  is  also  significantly  greater  than  Thorp's 
value. 

Under  the  assumption  that  the  anomaly  is  truly  a  relaxation  process,  the 
difference  between  the  two  relaxation  frequencies  can  be  attributed  to  the 
difference  in  temperatures.  However,  the  difference  in  coefficients  is 
somewhat  greater  than  the  estimated  experimental  error.  Despite  this 
disagreement,  a  second  relaxation  still  seemed  to  be  the  most  promising 
explanation  of  the  anomaly. 
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Figure  7:  Red  See  data  end  2-component  model. 

Browning  el.  el.  carried  out  a  number  of  sound-channel  experiments  ‘  j 
test  the  relaxation  hypothesis.  In  the  Red  Sea  experiment  of  1969  1 12],  the 
track  was  300  km  long  and  the  water  depth  was  750  m  or  greater.  Sound- 
velocity  profiles  showed  a  strong  sound  channel  with  axis  near  200  m.  The 
ambient  temperature  was  22®C  on  the  axis  and  30°C  at  the  surface. 
Measurements  were  made  exis-to-axls  using  SUS  charges  at  200  m  depth. 

Relaxation  frequency  can  be  expected  to  Increase  exponentially  with 
temperature  over  the  limited  range  of  temperatures  encountered  In  the 
sea.  Based  on  extrapolation  of  the  Thorp  and  Leroy  values,  the  expected 
relaxation  frequency  in  the  Red  Sea  experiment  was  approximately  2  kHz. 
However,  the  curve  fit  of  Figure  7,  Indicated  a  significantly  lower  value. 

The  Red  Sea  formula  used  In  the  data-flt  is  given  by: 

A=  0.008  f2  ♦  0. 1 48  fr  f2/(  fr2  ♦  f2 )  dB/km  (fr=  1 .5  kHz) 

The  anomaly  coefficient  is  clearly  in  better  agreement  with  Leroy  s  value 
than  that  of  Thorp  and  the  magnesium  sulfate  coefficient  is  significantly 
greater  that  predicted  by  the  S&M  formula. 

The  relaxation  hypothesis  remained  promising  despite  inconsistencies; 
however,  it  was  not  until  1973  that  it  was  firmly  established  by  an 
entirely  different  experimental  technique,  i.e.  the  T-jump  method. 
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T-Jump  Measurements 


Figure  8:  Temperature- jump  apparatus. 

In  1973  Yeager  et.  al.  1131  reported  results  of  T-jump  experiments  with 
both  synthetic  and  actual  sea-water.  In  this  method,  a  very  rapid  rise  in 
temperature  is  produced  by  discharge  of  the  capacitor  through  electrodes 
in  the  cell  containing  the  solution,  as  shown  in  Figure  8.  A  chemical  pH 
indicator  was  used.  Light  transmission  to  the  photocell  measures  the  pH 
change  produced  by  the  temperature  jump  and  the  relaxation  frequency  is 
determined  from  the  rate  of  change  as  measured  by  the  oscilloscope. 

Sea  water  showed  the  1  kHz  relaxation,  proving  that  the  equilibrium  is 
the  acid/base  type.  Synthetic  sea  water  (SSW)  was  tested  by  addition  of 
chemicals  to  distilled  water  according  to  the  formula: 

Constituent  Concentration  (mM) 

Sodium  Chloride  (NaCl)  400 

Magnesium  Sulfate  (MgS04)  35 

Magnesium  Chloride  (MgCl2)  20 

Calcium  Chloride  (CaCl2)  10 

Sodium  Bicarbonate  (NaHC03)  2.5 

Boric  Acid  (B(0H)3)  0.5 

where  concentrations  are  in  millimoles  per  liter  (mM).  The  constituent 
responsible  was  Identified  as  boric  acid.  The  relaxation  rate  was  found  to 
Increase  with  concentration  and  temperature  but  pH  had  negligible  effect. 

In  aqueous  boric  acid,  the  ionization  mechanism  is  believed  to  be. 

B(0H)3  ♦  OH'  e  B(0H)4" 

where  0H‘  is  the  hydroxyl  ion  and  B(0H)4-  is  the  borate  ion.  However,  the 
concentration  of  the  carbonate  ion  is  much  greater  in  sea  water  and  an 
exchange  reaction  was  proposed  as  the  more  likely  mechanism. 
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Figure  9:  B(0H)3  and  C02  specie  concentrations  vs  pH. 


The  proposed  Ionization  equilibrium  in  sea  water  is  given  by: 

B(0H)j  ♦  CO32-  ♦  HjO  B(0H)4"  ♦  HC0j“ 

For  chemical  purposes,  the  individual  equilibria  C8n  be  written: 

HCOj-  **  H4  ♦  CO/'  K^l  H4]|  COj^I/IHCOj-] 

B(0H)3  ♦  F^O  O  H4  ♦  B(0H)4-  Kg=l  H+]|B(0H)4-]/lB(0H)3J 

where  I  H4]  is  the  hydrogen  ion  concentration  and  pHs-log(I  H4)).  The 
equilibrium  constants  in  these  terms  are  pKBa9.2  and  pK^.B. 

Concentration  of  the  carbonate  Ion,  for  example.  Is  given  by  the  formula: 

I  CO^- 1 »  (C02l  1 0x/(  1+  I0X  )  x=pH-pKc 
where  [C021  is  the  total  carbon  dioxide  concentration. 

Possible  effects  of  ion-pairing  of  carbonate  with  metallic  ions  were 
also  pointed  out.  Carbonate  has  a  strong  affinity  with  sodium,  magnesium 
and  calcium  Jons.  Bicarbonate  affinities  are  much  smaller  and  can  be 
Ignored.  The  principal  equilibria  to  be  considered  are  therefore: 

Na4  ♦  C032’  «■>  NaCOj"  Ks  1.6/M 

Mg24  ♦  COj2' «  ttgC03  K«53/M 
Ca24  ♦  C032'  o  CaC03  K *  1 25/M 

In  sea  water,  most  of  the  carbonate  exists  In  essoci&s-'d  form  and  the 
C03  curve  of  Figure  9  represents  the  total  of  the  four  carbonate  species. 

Borate  ion  formulae  are  similar  but  the  affinities  are  much  smaller. 
Boric  acid  and  bicarbonate  affinities  are  negligible 


Seo-Doto  Analusis 


Temperature  *C 

Figure  10:  8(0H)5  relaxation  frequency  data. 

In  order  to  estimate  the  relaxation  frequency  of  boric  acid  with  the 
two-component  model,  the  measured  absorption  spectra  had  to  be  fitted  by 
a  curve  with  two  other  parameters;  namely,  the  magnesium  sulfate  and 
boric  acid  coefficients.  Because  of  limited  accuracy  of  the  low-frequency 
data,  the  boric  acid  relaxation  frequency  often  tends  to  be  difficult  to 
determine  with  the  required  accuracy. 

Figure  10  includes  only  those  experiments  for  which  reasonably  good 
estimates  of  relaxation  frequency  could  be  obtained  from  analysis  of  the 
experimental  data.  The  Atlantic  and  Red  Sea  experiments  were  believed  to 
be  the  most  accurate  and  were  heavily  weighted.  The  solid  line  is  given  by 
fr=0.9x10T/,o°(kH2),  which  intersects  both  points.  In  the  analysis  of  the 

other  data,  this  formula  was  taken  to  be  the  best  estimate  of  relaxation 
frequency  vs  temperature  at  that  time. 

Except  for  Leroy’s  Mediterranean  point,  there  is  reasonable  agreement 
within  the  limits  of  estimated  experimental  error.  Further  investigations 
have  resolved  all  of  the  dlscrepencies,  as  will  be  seen  in  the  later 
sections  of  this  paper. 
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Figure  1 1:  Boric  acid  coefficient  a(B)/YThorp  vs  pH. 

For  the  exchange  mechanism  proposed  by  Yeager  et.  al ,  the  boric  acid 
coefficient  a(B)  is  proportional  to  the  product  of  the  two  reactants  on 
either  side  of  the  equilibrium  as  Indicated  earlier.  Taking  boric  acid  and 
carbonate  as  the  reactants  we  have: 

|B(0H)31IC0j2']  »  IBHC02)  10*/(1  +  10*)(I*I0¥  ) 

where  iBl  is  total  boric  acid  concentration  and  x=pH-pKc,  y=pH-pKg 

The  coefficient  should  therefore  increase  exponentially  and  then  decay 
after  reaching  a  maximum  value  near  pHs9.5.  In  the  curve  Figure  1 I,  the 
coefficient  is  expressed  relative  to  Thorp's  value.  In  the  sea-water  pH 
range,  the  coefficient  Increases  approximately  as  lOpM.  In  practice,  pH>8.5 
Is  not  possible  In  sea  water  because  calcium  carbonate  will  precipitate 
out  of  solution. 

The  pH  dependence  of  the  experimental  sea-data  was  investigated  I  Ml. 
Analysis  was  based  on  the  assumption  that  the  rate  constant,  which 
controls  the  relaxation  frequency.  Is  the  main  temperature-dependent 
factor.  The  relaxation  amplitude  parameter,  e(B),  is  easily  determined 
when  the  relaxation  frequency  is  calculated  from  the  data-f it  of  Figure  10. 
The  value  should  also  be  nearly  independent  of  temperature  because  the 
rate  constant  then  cancels  out  of  the  equation. 


-  13- 


/ 


a(BX 


Figure  12:  B(0H)3  parameter  a(B)//Thorp  vs  pH 

All  sea  data  were  fitted  by  the  two-relaxation  model  l  Ml.  Values  for 
the  magnesium  sulfate  component  were  obtained  by  matching  the  spectra 
at  the  higher  frequencies.  Other  components  of  low-frequency  loss  were 
also  added  when  indicated,  (see  the  later  sections  on  diffraction  and 
scattering  loss).  The  parameters  of  the  boric  acid  spectrum  were  then 
determined  by  adjusting  the  amplitude  and  relaxation-frequency 
parameters  for ’best-fit’. 

When  the  relaxation  frequency  could  not  be  estimated  with  sufficient 
accuracy,  the  boric  acid  coefficient  Am*x=a(B)  fr  was  measured  The  value 

of  a(B)  was  then  obtained  by  dividing  Am«x  by  the  appropriate  relaxation 
frequency  from  the  data-fit  of  Figure  10.  The  values  shown  in  Figure  12 
are  relative  to  Thorp's  value  a(B)=0.l  1,  which  is  still  believed  to  be  the 
best  approximation  for  the  two-relaxation  model. 

The  standard  deviation  of  the  data  vs  pH  shows  a  correlation  of  better 
than  90%  and  no  temperature  trend  Is  evident,  supporting  the  assumption 
that  relaxation  frequency  is  the  only  temperature-dependent  parameter. 
Using  Thorp's  pH=8  as  the  standard,  the  approximation  for  the  absorption 
parameter  becomes  a(B)=0. 1  lxlO(pH‘9) 
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Resonator  Measurements 


Figure  13:  Resonator  apparatus  (ref.  16) 


The  temperature-jump  measurements  showed  that  only  boric  acid  can 
produce  a  relaxation  in  the  1  kHz  range.  However,  since  the  method  permits 
measurement  of  relaxation  frequency  and  not  absorption,  it  was  not  yet 
clear  what  part  it  played  in  the  actual  absorption  mechanism. 

Resonator  measurements  in  the  1  kHz  range  are  not  feasible  because  of 
size  limitations.  Subsequent  investigations  by  Simmons  1151,  using  the 
resonator  method  at  higher  frequencies  and  the  temperature-jump  method, 
proved  that  boric  acid  is  essential  to  the  low-frequency  sound-absorption 
process.  However,  since  aqueous  boric  acid  showed  comparable  relaxation 
frequency  as  well  as  absorption,  the  mechanism  remained  uncertain. 

From  1979-1983,  resonator  measurements  of  synthetic  sea  water  (SSW) 
were  carried  out  at  NUSC  by  Mellen  et.  al.  1161.  The  apparatus  is  shown  in 
Figure  13.  The  useable  frequency  band  was  limited  to  10-100  kHz  and  the 
errors  were  estimated  to  increase  from  *0.2  to  *1  dB/km  in  this  range. 

In  order  to  measure  the  relaxation  parameters,  the  reactions  had  to  be 
maneuvered  into  the  measurement  window  by  increasing  concentrations. 
Sea-water  relaxation  parameters  were  obtained  by  extrapolating  back  to 
normal  concentrations.  Sulfate  was  omitted  in  order  to  eliminate  the 
magnesium  sulfate  relaxation  and  thereby  improve  the  accuracy.  Other 
relaxations  involving  magnesium  were  also  discovered,  however,  of  these, 
only  the  magnesium  carbonate  relaxation  plays  any  significant  role  in  sea 
water  (161.  All  relaxations  were  measured  separately  and  in  combination 
to  determine  interaction  effects  between  them. 
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Figure  14:  Boric  ocid  absorption  spectra  at  pH=9. 


Magnesium  chloride  has  negligible  effect  on  the  boric  acid  relaxation 
and  was  omitted  in  the  experiments  of  Figure  14  in  order  to  eliminate  the 
undesired  relaxations.  In  the  top  spectra,  calcium  chloride  was  included  at 
10  mM  and  the  sodium  bicarbonate  concentration  was  reduced  to  1  mM  in 
order  to  prevent  precipitation  of  calcium  carbonate.  Relaxation  frequency 
inert  ,es  linearly  with  boric  acid  concentration  to  nearly  35  kHz.  The 
lower  spectra  (calcium  chloride  omitted)  show  the  same  effects  when  the 
sodium  bicarbonate  concentration  is  increased  end  the  concentration  of 
boric  acid  is  held  constant  (2  mM),  confirming  the  exchange  mechanism. 


I  10  100 

Concentration  Factor  (mn)*2 


Figure  15:  Boric  acid  Amax  det8  vs  concentration. 

in  figure  15,  the  maximum  values  of  each  spectrum  of  Figure  14  are 
plotted  vs  the  relative  concentration  of  the  varied  constituent. 

Am«  Is  proportional  to  the  product  of  the  reactants.  The  concentration 
factors  (mM)2  for  the  two  cases  are  given  by: 
lB(OH)j]  X  ICOjtotlconst. 

ICOj2’]  x  (B(0H),)cowt. 

WKh  pH«9  held  constant,  the  concentration  of  the  varied  reactant  is 
simp'y  a  constant  factor  times  its  total  concentration.  The  calcium 
carbonate  ion-pair  concentration  is  included  in  the  first  case.  Calcium 
chloride  was  omitted  In  the  second  case. 

Both  data  sets  are  proportional  to  the  concentration  factor;  however,  the 
absorption  increases  dramatically  when  cslcium  Is  included,  which  clearly 
Indicates  a  large  contribution  to  AV.  The  other  metallic  ions  do  not  show 
any  significant  effects. 

Extrapolated  to  sea-water  concentrations  and  pH,  the  data  with  calcium 
are  in  good  agreement  with  sea  measurements.  Without  calcium,  the  values 
are  too  low  by  almost  an  order  of  magnitude. 
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Figure  16:  Comparison  of  resonator  and  sea  data  vs  pH. 


Figure  16  shows  a  comparison  of  the  boric  acid  resonator  data  with  the 
sea  data  of  Figure  12.  The  resonator  data  were  obtained  at  elevated  boric 
acid  concentrations  and  wer  e  corrected  to  sea-water  concentration  by 
extrapolation. 

The  resonator  experiments  confirm.  In  detail,  that  the  boric/carbonate 
exchange  Is  the  dominant  mechanism  and  that  all  assumptions  are  valid. 

As  the  result  of  this  work,  it  Is  believed  that  the  boric  8cid  absorption 
mechanism  is  now  understood  to  almost  the  same  degree  cs  that  of 
magnesium  sulfate.  However,  to  construct  an  practical  absorption  model 
based  solely  on  relaxation  chemistry  would  be  far  too  complex  and  the 
accuracy  would  be  limited  to  that  of  the  supporting  measurements  as  well. 
Therefore,  the  combined  results  of  resonator,  T-Jump  and  sea  experiments 
have  been  taken  Into  consideration  In  the  development  of  a  simplified 
model  for  the  boric-acid  mechanism. 

Tnere  are  jO  number  of  possible  models  for  the  boric  acld/carbonate 
equIHbrlumjln  sea  water.  One  that  Is  both  tractable  and  consistent  with 
the  experimental  data  is  preferred  as  a  matter  of  convenience  and  may  not 
represent  all  the  essential  chemical  details  Involved. 
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lorlc  Acid  Mechanism 


B(0H)3+  CoC03+H20  B(OH);+  HCOg+Ca2+ 

principal  exchange  reaction 

B(OH)  ♦C03”+Ca2+H20  «-»  B(OH)^  *HCO~ +Ca2+ 

I3  2  I«  , 

B(0H)3+CaC03*H20  <- -»  CaB(0H)4  *HC03 
coupled  exchange  reaction 


frS2ff(IB(0H)3j  +  [HC°’1/K,)  (HZ) 


Am«X  ft 


(AV)2k>[B(0H)3HC0M.tlPC2(WD/-c,.) 

2R(T<-273) 

absorption  formula  parameters 


AV  s*-22ml/mole 

k>  s*4xl  06+T/70/mol*sec 
R=gas  constant  T=temp  *C 
P  =density  c=sound  speed 

_ 

Figure  17:  Boric  acid  relaxation  equilibria. 


The  top  box  of  Figure  17  shows  the  principal  equilibrium  describing  the 
essential  chemical  factors  involved.  However,  the  relaxation  frequency  for 
this  model  depends  on  calcium  concentration,  which  Is  not  observed 
experimentally. 

The  second  box  shows  details  of  the  proposed  ion-pair  equilibria,  in  this 
case,  the  two  sides  are  separated  by  equally  slow  steps  and  relaxation 
frequency  Is  quite  independent  of  calcium  concentration.  This  model  is 
used  as  the  basis  of  the  present  analysis. 

Simplified  formulae  for  the  absorption  parameters  for  the  model  are 
given  In  the  third  box.  Species  having  negligible  effect  on  relaxation 
frequency  are  omitted  In  the  formula. 

The  bottom  box  shows  the  estimated  values  of  the  relaxation  parameters 
AV,  k>  and  K,. 
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Other  Relaxations 


Frequency  (khte) 

Figure  18.  Absorption  spectra  of  B(0H)3/C02/Mg/Ca  system. 

The  contributions  of  two  other  relaxations  involving  magnesium  are 
shown  in  Figure  18.  The  mechanisms  of  both  appear  to  be  similar  to 
magnesium  sulfate  and  are  modeled  by  two-step  equilibria: 

(MgCOj),  o  (figC03)„  **  ng2+  ♦  C032' 

(MgB(0H)4+),  f*  (MgB(0H)4+)j ,  o  Mg2*  ♦  B(0H)4’ 

Both  relaxations  have  been  studied  independently.  In  the  magnesium 
carbonate  experiments,  boric  acid  was  omitted.  In  the  magnesium  borate 
experiments,  sodium  bicarbonate  W8S  omitted.  Without  carbonate,  boric 
acid  absorption  is  so  small  that  it  can  be  ignored. 

Temperature,  concentration  and  pH  dependencies  of  both  relaxations 
were  determined  from  the  relaxation  spectra. 

With  the  parameters  of  the  two  relaxations  known,  the  two  systems 
were  then  combined  to  determine  any  interaction  effects.  The  results 
shown  in  Figure  18  indicate  that  interactions  are  negligible  and  that  all 
three  components  are  simply  additive. 

For  normal  sea-water  boric  acid  concentration,  the  magnesium  borate 
component  becomes  negligible.  The  magnesium  carbonate  component,  while 
small,  can  become  significant  at  high  pH  values  [171. 
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Absorption  Model 


A=A1(MgS04)+A2(B(0H)3)+A3(MgC03) 

A„=  anf2f„/(f2+fn2) 

a,  =  0.5  xlo‘d(l!,7',/2a  f,=50*IOT/<° 

82=o..x,o<'h-8)  f2=oW'70 

8  3=  0.03x!0^7H"8^  fj=4.5xloY30 


Atlantic  4*C  pH  8.0 

A=0.007f  2  *0.1  f 2/(1»  f  2)+0. 18  f 2/(62>f 2) 
H.Pacific  4*C  pH  7.7 

A=0.007f2*0.05f2/(l+f2)+0  09  f2/(62+f2) 

Mediterranean  14*C  pH 8.3 

A=  Q.QQfif  2+0.26f  2/(1.42+f  2)+0.78  f2/(122+f  2) 

Red  Sea  22*C  pH  8.2 

A=0.0Q4f  2*0.27f  2/(1.82*f  2)  +  1.  If  2/(242 +f  2) 

sub-Arctic -1*C  pH  8.3 

A=0.01f2  +  0.17  f2/(0  852*f2)  +  0.24f2/(42*f2) 

Figure  19:  Simplified  absorption  formulae. 


The  simplified  formula  for  the  3-relaxation  model  is  shown  at  the  top  of 
Figure  19  where  A  is  the  absorption  coefficient  (dB/km),  f  is  frequency 
(kHz),  T  is  temperature  (°C)  and  pH=8.0  is  the  reference  value.  The  pure- 
water  term  is  omitted,  making  the  formula  valid  for  frequencies  less  than 
roughly  100  kHz.  The  pressure  factor  in  the  magnesium  sulfate  formula  of 
Fisher  and  Simmons  (181  for  has  been  modified  as  the  depth  term  D(km)  in 
this  model  Depth  dependencies  of  the  other  relaxations  are  not  yet  known, 
however,  measurements  in  both  deep  and  shallow  channels  indicate  that 
the  effect  on  boric  acid  must  be  very  small.  Magnesium  carbonate  effects 
may  be  greater  but  can  be  neglected  because  its  contribution  is  so  small. 
Salinity  dependence  is  neglected  temporarily  but  will  be  considered  in  a 
later  section. 

Specific  values  for  several  experimental  locations  are  also  shown  in  the 
bottom  box.  Note  that  the  magnesium  sul f ste  terms  A1  are  approximations 
valid  only  for  f«f,. 
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Model  and  Data 


Figure  20.  Thorp's  data  end  model. 

Figure  20  shows  Thorp's  data  compared  to  the  3-component  relaxation 
model.  Note  that  the  magnesium  sulfate  (A  I)  coefficient  is  0.007  rather 

than  the  earlier  value  0.009  used  in  the  2-component  model.  The  new  value 
is  based  on  a  very  careful  laboratory  re-evaluation  of  magnesium  sulfate 
absorption  in  sea  water  by  Fisher  end  Simmons  1 1 81. 

Note  also,  that  the  boric  acid  (A2)  coefficient  is  0.1  compared  to  0.1 1  in 
the  2-component  model.  The  magnesium  carbonate  (A3)  component  makes 
up  for  both  differences  and  the  fit  to  the  data  is  somewhat  belter. 

The  possible  existence  of  a  third  relaxation  in  sea  water  was  proposed 
earlier  by  Fisher  (19J.  He  observed  a  persistent  bump  in  the  data  around 
3-4  kHz  that  could  not  be  removed  by  any  consistent  readjustment  of  the 
parameters  of  the  2-component  model.  The  discovert;  of  the  magnesium 
carbonate  relaxation  in  Lake  Tanganyika  and  subsequent  study  of  the 
mechanism  put  the  hypothesis  on  very  frm  ground,  in  fact,  resonator 
measurements  of  sea  water  at  *.Hs8.5  revealed,  beyond  all  doubt,  the 
deficiencies  of  the  2-component  model.  Regardless  of  any  reasonable 
readjustment  of  the  parameters,  measured  values  were  much  too  high 
without  the  third  component  1171. 
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In  the  North  Pacific  model,  the  lower  value  pHs7.7  reduces  both  the  boric 
acid  (A2)  and  the  magnesium  carbonate  (A3)  coefficients  by  a  factor  of 
two  compared  to  Thorp.  Relaxation  frequency  depends  only  on  temperature 
and  remains  the  same.  Note  that  the  dB/km  scale  has  been  reduced  by  the 
factor  10  In  Figure  21. 

Sea  water  is  highly  alkaline.  Values  of  pH  In  the  World  Ocean  range  from 
about  7.7  to  8.3  Whereas  the  neutral  value  Is  7.  Surface  values  average 
around  8.1  and  tend  to  have  a  much  narrower  range  Evidently,  circulation 
and  contact  with  the  atmosphere  are  Involved  In  controlling  pH.  The  lowest 
values  are  found  In  deep  waters  where  the  circulation  times  are  the 
longest,  which  would  account  for  the  very  low  absorption  observed  in  the 
North  Pacific  sound-channel. 

At  the  higher  latitudes,  the  sound-channel  axis  comes  to  the  surface  as 
the  thermocllne  diminishes  and  the  absorption  increases  because  the  pH  is 
higher.  The  very  low-absorption  phenomenon  Is  therefore  probably  limited 
only  to  the  central  North  Pacific  Ocean.  I 
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Figure  22:  Red  Sea  data  and  model. 

The  Red  Sea  data  were  originally  fitted  with  a  2-relaxation  model  and 
the  results  Indicated  a  low-frequency  relaxation  at  1.5  kHz  compared  to 
the  formula  value  1.8  kHz.  Figure  22  shows  that  an  even  belter  fit  is 
achieved  with  the  3-relaxatlon  model  using  the  latter  value. 

The  major  remaining  anomaly  appears  to  be  the  1.7  kHz  relaxation 
frequency  of  Leroy's  Mediterranean  experiment.  Since  the  method  of 
measurement  is  different,  other  sources  of  error  must  be  considered. 

Experiments  reported  here  made  use  of  explosive  sources.  A  possible 
source  of  error  is  nonlinearity  wherein  the  overtaking  effect  regenerates 
some  of  the  lost  high-frequency  energy  at  the  expense  of  low.  Earlier,  an 
attempt  was  made  to  explain  the  entire  attenuation  anomaly  by  such  a 
mechanism.  However,  propagation  in  a  sound  channel  Involves  phase  shifts 
at  caustics,  which  reverses  nonlinear  effects  so  that  there  is  negligible 
accumulation  at  long  ranges.  Leroy's  ranges  were  limited  to  the  first 
convergence;  l.e  s40  km.  If  the  wavefront  remains  coherent  over  the  entire 
propagation  path,  effects  can  be  significant. 

The  recent  Mediterranean  sound-channel  experiment  appears  to  support 
this  argument  because  the  range  extended  to  more  than  600  km  and  results 
ere  consistent  with  the  model. 
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Figure  23:  Mediterranean  data  and  model. 

The  Mediterranean  experiment  was  carried  out  in  the  deep  Ionian  basin 
between  Sicily  and  Crete  In  July  1982  by  SACLANTCEN  111.  sound  velocity 
profiles  Indicated  a  strong  sound  channel  with  a  rather  broad  minimum  at 
100-200  m  depth.  Temperature,  salinity  and  sound-speed  profiles  were 
uniform  over  the  640  km  path.  Effects  of  the  rising  bottom  were  noted 
only  at  extreme  ranges  and  these  points  were  omitted  In  the  analysis. 

Measurements  of  pH  and  boric  acid  concentrations  also  showed  little 
variability  with  range  or  depth  over  the  track.  The  high  pH  value  is  typical 
of  surface  values  In  other  regions.  Both  the  boric  acid  and  carbon  dioxide 
concentrations  were  found  to  be  normal.  The  salinity  (S»38)  is  slightly 
higher  than  normal. 

At  pH=8.3,  the  boric  acid  coefficient  Is  a  factor  of  2  greater  than  that 
for  the  North  Atlantic.  However,  the  much  higher  temperature  Increases 
the  relaxation  frequency  so  that  the  absorption  spectrum  Is  not  increased 
by  such  a  large  factor  at  the  lower  frequencies. 

The  3-component  relaxations  calculated  for  the  experimental  conditions 
are  shown  in  Figure  23  and  the  fit  to  the  data  Is  clearly  excellent. 
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Figure  24:  Bismark  Sea  surface-duct  data  and  model. 

An  experiment  in  the  Bismark  Sea  by  Wylie  (201  indicates  the  need  for 
including  diffraction  loss  in  surface-duct  calculations.  In  this  case,  the 
duct  depth  was  approximately  90m.  Diffraction  effects  become  significant 
when  the  acoustic  wavelength  becomes  ito  large  for  total  refractive 
trapping.  The  consequence  is  an  additional  component  of  attenuation  that 
increases  rapidly  with  decreasing  frequency. 

Acoustic  trapping  is  the  result  of  an  isothermal  mixed-layer  at  the 
surface  overlying  the  thermocline.  The  sound-speed  gradient  in  the  layer 
can  be  approximated  as  that  in  deeper  isothermal  water.  The  diffraction 
loss  in  Figure  24  was  calculated  for  a  bilinear  gradient  with  the  normal 
positive  gradient  dcwn  to  the  duct  depth  and  a  negative  gradient  of  larger 
magnitude  below.  Note  that  trapping  is  never  total;  however,  the  loss 
decreases  exponentially  with  increasing  frequency  and  rapidly  becomes 
negligible.  The  only  effect  of  diffraction,  in  this  case,  appears  to  be  the 
absence  of  data  at  frequencies  below  effective  cutoff. 

The  absorption  model  was  calculated  for  the  ambient  temperature  3C''C 
and  the  pH=8.3  was  estimated  from  archival  data.  Agreement  with  the  data 
above  1  kHz  indicates  that  surface  scattering  was  negligible  as  expected 
for  the  low  sea  state 


Figure  25.  Gulf  of  Aden  surface-duct  data  and  model. 

In  the  Gulf  of  Aden  experiment  112],  the  duct  depth  was  only  30m  end  the 
calculated  diffraction  effect  remains  strong  to  well  above  1  kHz.  However, 
the  data  shows  increasing  attenuation  over  only  a  limited  range  and  the 
trend  reverses  as  frequency  decreases.  This  is  the  result  of  the  leaky 
surface-duct  signals  returning  from  the  deep  channel. 

The  absorption  model  was  calculated  for  the  ambient  temperature  30°C 
and  the  pH  value  estimated  from  archival  data.  Agreement  with  the  data 
above  1  kHz  indicates  plausibility  of  both  the  absorption  and  diffraction 
models. 

See  states  in  both  experiments  were  very  low  and  the  absence  of  surface 
loss  was  not  unexpected.  In  most  sound-channel  experiments,  refractive 
trapping  by  the  thermocline  tends  to  make  effects  of  surface  scattering 
negligible.  However,  if  the  thermocline  is  weak  or  absent,  the  excess  loss 
should  certainly  become  appreciable  at  moderate  sea  states.  In  such  cases, 
a  marked  frequency-dependence  is  expected,  particularly  at  the  lower 
frequencies.  For  acoustic  wavelengths  large  compared  to  the  surface 
wavelength,  the  excess  loss  should  fall  off  very  rapidly. 
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Lake  Experiments 


The  experimental  program  of  Browning  et.  a).  1121  included  two  lake 
experiments  for  comparison  purposes.  Both  gave  surprisingly  anomalous 
results  compared  to  pure  water  as  seen  in  Figure  26. 

Attenuation  in  Lake  Tanganyika  was  found  comparable  to  the  magnesium 
sulfate  relaxation  in  sea  water.  Analysis  of  the  constituents  showed  that 
the  concentration  was  too  low  by  a  factor  of  30  to  account  for  the  result; 
however,  both  the  carbon  dioxide  content  and  pHz8.5  were  abnormally  high. 

/  Chemical  synthesis  of  artificial  lake  water  in  other  resonator  experiments 

showed  the  mechanism  to  be  the  magnesium  carbonate  relaxation. 

Attenuation  in  Lake  Superior  of  W8S  found  to  be  comparable  to  that  of 
boric  acid  in  sea  water.  In  this  case,  the  anomaly  has  been  attributed  to 
internal  scattering.  Frequency-independent  attenuation  is  apparently  the 
result  of  multiple-scattering  by  large-scale  temperature  inhomogeneities 
within  the  thermocline.  Energy  gradually  diffuses  out  of  the  channel  by  a 
“random-walk"  process  and  is  absorbed  in  the  bottom.  This  mechanism 
seems  to  be  common  to  all  sound-channel  experiments;  however,  it  is 
.  „  usually  observed  only  at  very  low  frequencies  (211.  Much  of  the  spread  in 

scattering  loss  can  be  attributed  to  differences  in  channel  strength; 
however,  no  reliable  method  of  prediction  has  yet  been  found. 
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Scattering  Effects 


Figure  27:  Hudson  Bay  data  and  model. 

Hudson  Bay  1 12]  was  the  first  sea  experiment  to  show  excess  frequency- 
independent  loss  and  served  to  put  the  scattering  hypothesis  on  more  firm 
footing.  The  puzzling  result  of  the  earlier  Lake  Superior  experiment,  in 
particular,  was  clearly  attributable  to  the  same  mechanism. 

Since  the  loss  is  independent  of  frequency,  ray-diffusion  theory  will  be 
applicable  and  losses  can  be  calculated  from  the  diffusion  equation  121). 
Diffusion  loss  depends  inversely  on  the  strength  of  the  channel,  in  order  to 
reach  the  bottom,  axial  rays  must  overcome  a  potential  barrier  AC,  which 
Is  the  difference  In  sound-speed  between  the  bottom  and  channel  axis.  For 
small  angles,  loss  is  Inversely  proportional  to  the  square  of  the  critical 
angle  of  the  bottom-grazing  ray. 

In  Figure  27,  the  dashed  line  is  the  estimated  total  absorption,  end  the 
excess  scattering-loss  is  taken  as  0.04  dB/km  independent  of  frequency. 
Water  depth  in  Hudson  Bay  is  only  100m  and  trapping  is  weak.  Although  the 
Lake  Superior  depth  is  nearly  the  same,  a  slight  negative  temperature- 
gradient  reduces  the  channel  strength,  making  the  loss  approximately  four 
times  greater. 
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Figure  28:  Gulf  of  Aden  data  and  model. 

In  the  Gulf  of  Aden  experiment  (12),  there  was  a  strong  Intermediate 
sound  channel  caused  by  outflow  of  water  from  the  Red  Sea.  This  is  in 
addition  to  a  deeper  channel  and  the  surface  duct  discussed  earlier.  Source 
and  receiver  were  located  In  the  Intermediate  channel.  The  experimental 
track  was  500  km  in  length  and  the  total  water  depth  was  roughly  1  km. 
Absorption  parameters  are  estimated  for  the  ambient  temperature  13*0 
and  archival  value  pH=7.9  at  axis  depth. 

Despite  the  relatively  strong  trapping  in  the  Intermediate  channel,  there 
Is  marked  evidence  of  scattering  loss.  This  may  be  the  result  of  instability 
in  the  upper  part  of  the  thermocline  where  the  temperature  falls  from 
30*C  at  the  surface  to  20’C  at  100  m.  Sound-speed  profiles  Indicated  a 
great  deal  of  fluctuating  fine-structure  and  this  could  account  for  the 
reverse  scattering  back  Into  the  surface  duct  as  noted  in  Figure  25.  In 
addition,  scattering  into  the  deeper  channel  may  also  be  a  significant 
factor. 

In  Figure  28,  an  additional  component  of  0.02  dB/km  has  been  added  to 
fit  the  data. 
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Figure  29:  Baltic  Sea  data  and  model. 

Interpretation  of  attenuation  measurements  In  the  Baltic  Sea  (221 
Indicates  need  for  estimating  sound  absorption  in  low-salinity  waters. 

While  the  magnesium  sulfate  component  depends  simply  on  the  first 
power  of  concentration,  the  low-frequency  boric  acid  relaxation  is  more 
complex.  The  mechanism  involves  exchange  with  carbonate  in  which  the 
Ion-pair  calcium  carbonate  governs  the  absorption  magnitude. 

Measurements  of  absorption  by  the  resonator  method  show  roughly  1/2 
power  dependence  on  calcium  concentration  In  the  range  of  Interest. 

The  Baltic  Sea  salinity  of  8  is  about  25*  normal  while  the  carbonate  and 
acidity  are  both  normal.  Since  loss  is  multiplicative,  S1 5  dependence  is 
Indicated.  This  would  yield  a  reduction  in  magnitude  by  the  factor  8. 
However,  a  more  exact  calculation,  taking  into  account  the  coupled-ion 
effects,  shows  somewhat  less  reduction.  In  addition,  relaxation  frequency 
is  lower.  T-Jump  measurements  of  Baltic  water  samples  by  Candau  1231 
(unpublished)  show  that  it  is  only  0.6  kHz,  which  agrees  with  theory.  All 
factors  have  been  Included  In  the  calculations. 

Comparison  of  the  measured  data  in  Figure  29  shows  that  attenuation  is 
still  anomalously  high  by  a  large  factor.  The  relative  weakness  of  the 
lower  part  of  the  channel  due  to  the  small  water-depth  (*  1 00  m),  is  again 
indicatatlve  of  loss  by  internal  scattering. 
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Figure  30:  Tasman  Sea  data  and  model. 

In  the  deep  ocean,  the  sound-channel  is  much  stronger  and  scattering 
loss  ranges  from  negligible  to  only  a  few  dB  per  thousand  kilometers.  The 
measurements  in  the  Tasman  Sea  by  Bannister  et.  al.  124)  in  Figure  30  are 
typical.  In  both  experiments  shown,  measurements  were  made  In  the  deep 
channel  and  the  maximum  range  was  roughly  3000  km.  The  absorption 
spectrum  was  calculated  for  pH=7.9,  a  value  estimated  from  archival  data. 

In  the  lower  curve,  a  scattering  loss  of  3x|l O'4  dB/km  has  been  added  in 
the  data-flt.  This  amounts  to  only  1  dB  total  over  the  path  and  Is  probably 
insignificant.  In  the  upper  data,  the  scattering  loss  is  taken  as  2xl0'3 
dB/km  and  amounts  to  roughly  6  dB  total,  which  is  considered  well  beyond 
the  limits  of  experimental  error. 

in  the  top  experiment,  the  track  encountered  the  Sub-Antarctic  Front 
where  the  channel  axis  rises  close  to  the  surface.  Range  dependence  of  the 
sound-speed  profile  may  account  for  some  of  the  anomaly;  however,  the 
absence  of  frequency  dependence  of  suggests  that  internal  scattering  is 
again  the  principal  mechanism. 

Regional  dependence  of  scattering  loss  has  also  been  investigated  in  the 
North  Pacific  by  Kibblewhlte  et.  al.  125).  Their  results  indicate  a  similar 
increase  at  higher  latitudes  that  the  range  dependence  of  the  sound-speed 
profile  does  not  account  for. 

i 
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Figure  31:  Sub-Arctic  data  and  model. 

Figure  31  shows  results  from  two  experiments  by  DREA  in  ice-free 
waters  west  of  Greenland.  Both  of  the  tracks  analyzed  were  in  excess  of 
300  km  in  length  with  water  depths  of  1-2  km.  The  sound-speed  profiles 
showed  a  distinct  channel  with  axis  depth  near  100  m  and  any  surface- 
scattering  effects  are  therefore  believed  to  be  negligible.  Absorption  has 
been  calculated  for  pH=8.2,  which  is  consistent  with  archival  data. 

The  southern  track  In  the  Davis  Strait  shows  clear  evidence  of  an  excess 
scattering-loss  component  despite  the  strong  channel.  The  data  are  fitted 
by  adding  the  value  8x10‘3  dB/km,  which  Is  comparable  to  the  highest 
values  reported  in  any  of  the  experiments  1251. 

The  northern  track  in  Baffin  Bay  also  shows  a  similar  excess  loss  but 
only  at  the  higher  frequencies.  At  100  Hz  It  becomes  negligible.  A  possible 
explanation  Is  that  the  bottom  absorption  at  lower  frequencies  may  be 
smaller  in  that  area,  causing  energy  to  be  reflected  back  into  the  channel. 
Alternatively,  the  Inhomogeneities  may  be  smeller  in  scsle  Scattering 
loss  can  be  expected  to  decay  rapidly  when  the  dimensions  become  small 
compared  to  the  acoustic  wavelength. 
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Figure  32:  Experimental  absorption  spectra  comparison. 

The  absorption  spectra  for  the  North  Atlantic,  North  Pacific  and  Eastern 
Mediterranean  experiments  are  compared  In  Figure  32.  The  coefficients 
have  been  divided  by  frequency  In  kHz  In  order  to  compress  the  scale. 

The  North  Atlantic  and  North  Pacific  are  apparently  the  extreme  cases 
for  deep  (>lkm)  waters  of  the  World  Ocean.  For  latitudes  less  than  45*,  the 
axis  depth  of  the  deep  sound-channel  is  roughly  1  km;  therefore,  the  two 
cases  represent  the  limits  of  variability  expectod  for  sound-channel 
propagation  In  this  latitude  range.  For  the  pH  range  7.7-0.O,  the  max/mln 
ratio  Increases  uniformly  with  decreasing  frequency,  approaching  a  factor 
of  two. 

The  Eastern  Mediterranean  is  the  highest  pH  so  far;  however,  since  the 
temperature  is  also  high,  the  curve  is  displaced  upward  in  frequency  and 
the  max/mln  ratio  is  slightly  less  than  a  factor  of  two  compared  to  the 
North  Atlantic. 

From  all  the  experimental  evidence,  it  appears  that  the  North  Pacific  and 
Eastern  Mediterranean  absorption  spectra  encompass  nearly  all  the  range 
of  expected  variability. 
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Figure  33.  Absorption  spectra  A/F  vs  pH  forT=4°C. 

The  principal  results  of  the  3-relaxation  model,  are  summarized  in 
Figures  33  and  34.  The  absorption  coefficient  A  has  again  been  divided  by 
the  frequency  in  kHz  in  order  to  compress  the  range.  The  predicted  spectra 
vs  pH  for  4®C  and  normal  salinity  (S=35)  are  shown  in  Figure  33  and  cover 
the  approximate  sea- water  range  of  pH  7.7-8  3  in  steps  of  O.l. 

The  effect  of  pH  on  the  absorption  spectrum  below  10  kHz  Is  mainly  due 
to  the  coefficient  of  the  boric  acid  component.  Below  1  kHz,  the  max/min 
ratio  approaches  the  factor  four. 

From  the  present  analysis  of  the  experimental  data,  the  overall  error  of 
the  absorption  model  Is  estimated  to  be  roughly  of  the  order  of  +15*.  In 
other  words,  the  RMS  error  of  an  estimated  coefficient  In  dB/km  at  any 
given  frequency  Is  not  expected  to  be  greater  than  this  amount  if  not 
limited  by  the  accuracy  of  the  environmental  factors.  Since  a  change  of 
0.05  pH  units  corresponds  to  a  )2%  change  in  the  boric  acid  coefficient, 
the  pH  errors  must  be  within  these  limits  in  order  to  realize  the  model 
accuracy.  The  major  limiting  factor  therefore  appears  to  be  the  accuracy 
and  availablility  of  of  the  pH  data-base. 
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Figure  34:  Absorption  spectra  A/F  vs  temperature  for  pH8.0. 

Figure  34  shows  the  predicted  spectral  dependence  on  temperature  for 
pH=8.0  and  normal  salinity  (S=35)  end  cover  the  nominal  sea-water  range 
0-30°C  in  5®C  steps. 

The  temperature  effect  in  the  model  comes  only  through  the  relaxation 
frequency  coefficients  of  all  three  terms.  Increasing  temperature  shifts 
the  curves  upward  in  frequency. 

Model  errors  involving  temperature  coefficient  are  also  estimated  to  be 
of  the  order  of  +15S.  Since  the  temperature  sensitivity  is  rather  low,  the 
value  need  only  be  known  only  to  within  roughly  ±2°C  to  realize  the  same 
estimated  accuracy. 

Salinity  corrections  have  been  omitted  In  the  simplified  model.  Errors 
are  probably  within  the  accuracy  limits  for  the  normal  sea-water  range 
32-38  ppl.  Sellnity  variations  are  evidently  the  result  of  evaporation  and 
dilution  by  fresh  water.  Therefore,  the  concentrations  of  all  constituents, 
except  carbon  dioxide,  tend  to  remain  in  constant  ratio.  A  multiplying 
factor  S/35  for  all  three  components  may  be  Justified  over  a  limited 
range,  say  30-40  ppt.  Beyond  this  range,  corrections  must  be  made  for 
changes  in  the  boric  acid  relaxation  frequency.  In  more  extreme  cases  like 
the  Baltic  See,  chemical  analysis  of  the  constituents  may  be  required. 
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Figure  35:  Absorption  ratio  A/Ao  vs  pH  and  temperature. 

In  Figure  35,  the  spectra  of  Figures  33  and  34  are  expressed  as  the  ratio 
A/Ao  where  the  Ao  is  calculated  for  the  North  Atlantic  parameter  values 
pH=8  0  end  T=4°C.  The  range  of  pH  is  7.7-8  3  in  steps  of  0. 1  and  the  range 
of  temperature  is  0-30°C  in  steps  of  5°C.  The  curves  indicate  the 
accuracies  required  to  realize  the  desired  error  limits. 
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Figure  37:  Pacific  Ocean  pH  profiles. 


Figure  37  shows  plots  of  typical  pH  profiles  for  a  range  of  latitudes  in 
the  Pacific  Ocean  from  GEOSECS  data  [261. 


\ 


Figure  38:  Indian  Ocean  pH  profiles. 


Figure  38  shows  plots  of  typical  pH  profiles  for  a  range  of  latitudes  In 
the  Indian  Ocean  from  GEOSECS  data  126). 

A  striking  feature  of  all  the  pH  profiles  is  the  relative  constancy  of  the 
values  at  the  surface  and  at  great  depths.  The  greatest  variability  occurs 
In  the  thermocline  region  around  1  km.  depth. 

For  surface-duct  propagation,  depth  variability  should  be  small  because 
both  the  pH  and  temperature  obviously  are  quite  constant  in  the  surface 
mixed-layer. 

For  convergence-zone  propagation,  most  of  the  travel  time  Is  In  the  deep 
water  where  the  pH  is  relatively  constant.  In  the  Southern  Ocean,  the  value 
Is  7.9.  In  the  North  Pacific,  pH  decreases  to  7.8.  In  the  North  Atlantic,  It 
increases  to  8.1  at  high  latitudes.  Variability  with  latitude  can  therefore 
be  expected  to  bs  relatively  easy  to  estimate. 

The  sound-channel  axis  depth  is  roughly  1  km  at  mid-latitudes  and 
rises  to  the  surface  at  high  latitudes.  The  net  result  is  that  the  pH  value 
on  the  axis  will  tend  to  Increase  at  the  higher  latitudes.  The  axial  pH  range 
is  approximately  7.7-8  3  and  values  tend  to  vary  rapidly  and  unpredlctably 
with  latitude.  Sound-channel  absorption  should  therefore  be  the  most 
difficult  to  estimate  with  acceptable  accuracy. 


Figure  39:  PH  test  profiles  and  ray  paths. 

In  order  to  estimate  the  effects  of  the  pH  profiles  at  mid-latitudes,  the 
two  typical  cases  illustrated  at  the  top  of  Figure  39  will  be  examined.  The 
sound-velocity  profile  (SVP)  is  typical  of  latitudes  less  than  45°.  Profile 
pH#  I  is  the  simple  exponential  decay  and  pH#2  is  the  case  where  there  is 
a  sharp  minimum  above  the  axis  as  in  the  North  Pacific.  Effective  values  of 
A  for  the  two  profiles  are  calculated  by  numerical  integration  of  the  loss 
over  single  cycles  of  the  sound-channel  and  convergence-zone  (CZ)  ray 
paths  shown  at  the  bottom. 
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Figure  40:  Absorption  ratio  A/ A*. 
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Figure  41:  High  latitude  test  profiles  and  ray  paths. 

In  order  to  estimate  effects  of  the  pH  profiles  at  higher  latitudes,  the 
two  cases  shown  at  the  top  of  Figure  41  will  be  examined.  The  SVP  is 
typical  of  latitudes  around  60°  Profiles  pH#3  and  pH#4  represent  current 
estimates  for  the  North  Atlantic  and  North  Pacific,  respectively.  The  ray 
paths  used  in  the  numerical  integration  of  loss  are  shown  below.  Note  the 
relative  lack  of  distinction  between  sound-channel  and  convergence  modes 
for  shallow  ray-angles.  At  higher  angles,  surface  reflection  is  involved. 
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Figure  42:  Absorption  ratio  A/A*. 

Figure  42  shows  the  effective  values  of  A/A*  where  A*is  calculated  for 
tjhe  2  km  depth  values  of  temperature  and  pH.  In  this  case,  errors  diminish 
fls  ray  angles  increase  while  the  opposite  occurs  if  axial  values  are  used. 

At  high  latitudes  in  the  Southern  Ocean,  profile  pH#  1  is  the  appropriate 
model  and  errors  for  this  case  are  roughly  doubled  compared  to  the  case  of 
profile  pH#3.  Interpolation  between  the  2  km  and  axial  values  in  the  ray 
analysis  would  clearly  help  to  keep  the  errors  within  acceptable  limits. 
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Global  Model 


A= A,  (MgS04  )+A2(B(0H)3)+A3(MgC03 ) 

An=(s/35)  a„f 2 

8j  =  0.5xl0-D<kn,)/20 
02s:0.1K 
03=  0.03  K 

Figure  43:  Global  model  absorption  formula. 

Variability  of  pH  is  clearly  the  major  limiting  factor  in  the  accuracy  of 
the  absorption  formula.  In  Figure  43,  the  pH  parameter  K=10(pH"8)  has  been 
substituted  in  the  simplified  absorption  formula.  Salinity  dependence  has 
been  taken  as  S/35  with  the  caveats  noted.  An  interim  model  is  proposed 
for  estimating  effective  K  values  for  specific  propagation  modes. 

Based  on  the  analysis  of  the  previous  section,  the  recommended  method 
is  to  use  axial  values  of  K  and  temperature  T  (°C)  for  the  sound-ch8nne] 
mode  and  the  2  km  depth  values  for  the  CZ  mode.  Errors  at  higher  latitudes 
con  be  reduced  by  interpolating  between  the  two  values,  depending  on  the 
ray  paths  involved.  Temperature  is  not  a  problem  since  it  can  be  easily  be 
derived  from  the  SVP  if  not  directly  available  from  XBT  data. 

Contours  of  pH  for  the  surface  and  for  the  0.5  km  and  1  km  depths  are 
shown  in  the  World  Ocean  Atlas  of  Gorshkov  1271.  Contours  for  2  km  depth 
are  also  included  in  Vol.  2.  Since  the  contour  intervals  are  0.1  pH  unit, 
interpolation  is  required.  Correction  to  in-situ  pressure  is  8lso  necessary. 

Pressure-corrected  pH  data  from  the  GEOSECS  report  have  also  been  used 
in  the  analysis,  allowing  estimation  of  the  2  km  contours  in  regions  of  the 
Pacific  Ocean  not  covered  by  the  Russian  work.  Discrepencies  between  the 
data  sets  in  other  regions  have  been  arbitrarily  resolved  by  adjustment  of 
the  contours  so  as  to  minimize  effects  on  estimation  error.  Considerable 
subjectivity  is  obviously  involved. 

The  K  contours  on  the  sound-channel  axis  in  Figure  44  are  based  on  an 
analysis  of  the  Russian  pH  contours  by  Lovett  (281.  The  2  km  depth  and 
surface  contours  derived  from  the  Russian  report  are  shown  in  Figures  45 
and  46,  respectively.  All  have  been  modified  in  taking  the  GEOSECS  dat8 
into  account. 
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Figure  44:  Sound-channe 
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Figure  46:  Surface  K  conlours.  j 
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Recommendations 


The  three-component  absorption  model  appears  to  be  firmly  established 
by  NUSC  laboratory  and  sea  experiments.  In  addition,  independent  studies 
at  the  Academia  Sineca  laboratory  confirm  the  results  in  detail  (29,30). 

The  proposed  simplified  model  has  been  shown  to  give  very  satisfactory 
field  predictions  within  the  error  limits  Imposed  bly  the  environmental 
data.  Future  refinements  may  require  slight  adjustment  of  the  parameters; 
however,  .ne  overall  accuracy  seems  to  be  adequate  for  current  purposes. 

The  main  limitation  of  the  global  model  is  the  accuracy  of  the  pH  values. 
Depth-dependence  does  not  appear  to  be  the  critical  factor.  The  simplified 
model  uses  axial  values  for  the  sound-channel  mode  and  2  km  depth  values 
for  the  other  modes.  Approximation  errors  are  evidently  smaller  than  the 
uncertainty  in  the  ofjjsolute  values  except  at  the  higher  latitudes  where 
interpolation  between  2  km  and  axial  values  will  help  improve  accuracy. 
The  K  contours  provided  in  this  report  should  give  reasonably  accurate 
approximations  for  absolute  pH  effects  in  most  of  the  World  Oceon.  The 
possible  exceptions  ere  regions  where  significant  discrepencies  between 
Russian  and  GEOSECS  values  have  been  noted.The  Russian  date  represent 
local  averages  of  archival  data  end  tend  to  be  lower  and  more  variable. 
Errors  involved  in  interpolating  the  data  can  also  be  significant. 

The  scope  of  this  report  has  been  limited  to  development  of  a  suitable 
absorption  model,  e  method  of  predicting  the  pH  effects,  and  a  rationale 
for  incorporation  intb  existing  propagation  models.  Analysis  has  shown 
that  the  estimated  overall  accuracy  ±15*  of  the  absorption  model  can  only 
be  realized  if  the  pH  error  is  ±0.05  units  or  less  end  this  rernoins  the 
crucial  factor.  Portions  of  the  GEOSECS  pH  data  have  been  used  to  analyze 
profile  effects,  to  Supplement  the  missing  Russian  2  km  contours  and  to 
modify  the  others  where  indicated.  Significant  disagreements  with  the 
Russian  data  have  been  noted.  It  is  therefore  recommended  that  the  entire 
GEOSECS  data  set  be  analyzed  for  the  purpose  of  either  refining  the  K 
contours  or  devising  a  more  accurate  data-base  format. 

At  sea,  the  simplest  solution  is  to  me-oure  surface  pH.  Measurement  of 
the  pH  profile  along  with  the  normal  XBT  would,  of  course,  help  to  resolve 
all  the  uncertainties.  Since  no  suitable  apparatus  seems  to  be  available  for 
this  purpose,  development  of  an  expendible  pH  prob*  may  be  indicated. 

Surface  scattering  is  another  loss  mechanism  that  can  be  important  in 
surface-duct  and  other  RSR  propagation  modes.  Previous  estimates  of 
scattering  loss  vs  weveheight  have  been  based  largely  on  surface-duct 
experiments.  Inaccurate  absorption  models  may  cause  the  effects  to  be 
seriously  overestimated  and  reassessment  using  the  improved  absorption 
model  is  also  recommended. 
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